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> 1 Introduction 

(N ' 

^ ! Hadronic jets are amongst the most striking phenomena in high-energy physics, and 

CN [ their importance is sure to persist as searching for new physics at hadron colliders 



becomes the main activity in our field. Signatures involving jets almost always have 
the largest cross sections but are the most difficult to interpret and to distinguish 
from background. Insight into the properties of jets is therefore doubly valuable: 
Oh! both test of our understanding of strong interaction dynamics and as a tool for 

Q-t, extracting new physics signals in multi-jet channels. 

In the present talk I shall concentrate on jet fragmentation and hadronization, the 
topic of jet production having been covered admirably elsewhere (0], |2|. The terms 
fragmentation and hadronization are often used interchangeably, but I shall interpret 
the former strictly as referring to inclusive hadron spectra, for which factorization 
'theorems'^ are available. These allow predictions to be made without any detailed 
assumptions concerning hadron formation. A brief review of the relevant theory is 
given in sect. |2|. 

Hadronization, on the other hand, will be taken here to refer specifically to the 
mechanism by which quarks and gluons produced in hard processes form the hadrons 
that are observed in the final state. This is an intrinsically non-perturbative process, 
for which we only have models at present. The main models are reviewed in sect. [3|. 
In sect. |] their predictions, together with other less mo del- dependent expectations, 
are compared with the latest data on single-particle yields and spectra. 

One of the most important objectives of jet studies is to understand the differences 
between jets initiated by different types of partons, especially quarks versus gluons, 

1 Plenary talk at XIX International Symposium on Lepton and Photon Interactions at High 
Energies, Stanford University, August 1999. 
2 Permanent address. 

3 I put this term in inverted commas because proofs of factorization in fragmentation do not really 
extend beyond perturbation theory H. 
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which could also be valuable in new physics searches. The wealth of recent data on 
quark-gluon jet differences is discussed in sect. |5|. Another goal is to compare jets 
produced in different processes such as e + e~ annihilation and deep inelastic scattering 
(DIS). This is being done by the HI and ZEUS collaborations at HERA; some of their 
results are discussed in sect. |6|. 

A good understanding of heavy quark jets is especially important since these 
are expected to be copiously produced in new processes such as Higgs boson decay. 
Section [7| discusses some recent results on heavy quark fragmentation. 

Finally, sects. § and || deal with correlation effects, first those of Bose-Einstein 
origin and then those expected in fully hadronic final states from e + e~ — > W + W~. 



A summary of the main points is given in sect. [10 



Regrettably I have no space to discuss many other relevant and interesting topics, 
such as: fragmentation function parametrizations |4]]; photon fragmentation function; 
polarization in hadronization; event shapes and power corrections ||; fluctuations 
and intermittency; dynamical two-particle correlations; heavy quark production in 
jets (g — > cc, bb); transverse energy flow in deep inelastic scattering; underlying event 
in DIS and hadron-hadron collisions; identified particle production in DIS; tests of 
QCD coherence; jet profiles and substructure |]. Apologies also to all those whose 
work I have omitted or mentioned only superficially. 

In citing the latest experimental data (mostly still preliminary) I have relied on 
the reference system of the EPS-HEP99 conference in Tampere, Finland ||, since 
few experimental papers were explicitly submitted to LP99. In the case of the large 
collaborations, papers submitted to EPS-HEP99 and/or LP99 can usually be found 
easily via the collaboration web pages 0H|~| . 



2 Jet fragmentation — theory 

Let us start by recalling the basic factorization structure of the single-particle inclusive 
distribution, e.g. in e + e" — > hX (fig. p]): 

F h (x, s) = J2 —Cfa a s (s))D?(x/z, s) 
i Jx z 

s = q 2 , x = 2p h ■ q/q 2 = 2E h / E cm 

where C\ are the coefficient functions for this particular process (including all selection 
cuts etc.) and D\ is the universal fragmentation function for parton % hadron h. 

The fragmentation functions are not perturbatively calculable but their s-dependence 
(scaling violation) is given by the DGLAP equation: 
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Figure V. Factorization structure of e + e — > ftJT. 



Thus they can be parametrized at some fixed scale So an d then predicted at other 
energies |@. 

In certain kinematic regions, higher-order corrections are enhanced by large log- 
arithms, which need to be resummed. At small x, log a; enhanced terms can be 
resummed by changing the DGLAP equation to 

This is commonly known as the modified leading-logarithmic approximation (MLLA) 
[|15| , |16| , [HJ . The effect of resummation is to generate a characteristic hump-backed 
shape in the variable £ = ln(l/x), with a peak at £ p ~ |lns. 

Large logarithms of ratios of invariants may also appear inside the coefficient 
functions Cj, for example in three-jet events when the angles between jets become 
small. In some cases these can be absorbed into a change of scale in the fragmentation 
functions. Examples will be encountered in sects. ^ and [|. 

Although universal, fragmentation functions are factorization scheme dependent. 
The splitting functions Pji are also scheme dependent in higher orders. To specify the 
scheme requires calculation of the coefficient functions to (at least) next-to-leading 
order. This has only been done in a few cases. Thus there is need for theoretical work 
to make full use of the data on fragmentation functions. 



3 Hadronization Models 



3.1 General ideas 



Local parton-hadron duality Hadronization is long-distance process, involving 

only small momentum transfers. Hence the flows of energy-momentum and flavour 
quantum numbers at hadron level should follow those at parton level. Results on 
inclusive spectra and multiplicities support this hypothesis. 

Universal low-scale a s fL9|, f20| , |21fl . Perturbation theory works well down to low 
scales, Q ~ 1 GeV. Assume therefore that a s (Q 2 ) can be defined non-perturbatively 
for all Q, and use it in evaluation of Feynman graphs. This approach gives a good 
description of heavy quark spectra and event shapes. 
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Figure 2: Cluster and string hadronization models. 



3.2 Specific models 



The above general ideas do not try to describe the mechanism of hadron formation. 
For this we must so far resort to models. The main current models are cluster and 
string hadronization. We describe briefly the versions used in the HERWIG and 
JETSET event generators, respectively. 



Cluster model p2|-p6 

qq, after the parton shower 
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The model starts by splitting gluons non-perturbatively, 
Colour- singlet qq combinations have lower 



masses and a universal spectrum due to the preconfinement |27, 2B property of 
the shower (fig. pD|). These colour-singlet combinations are assumed to form 



clusters, which mostly undergo simple isotropic decay into pairs of hadrons, cho- 



sen according to the density of states with appropriate quantum numbers [23 



This model has few parameters and a natural mechanism for generating trans- 
verse momenta and suppressing heavy particle production in hadronization. 
However, it has problems in dealing with the decay of very massive clusters, 
and in adequately suppressing baryon and heavy quark production. 



String model 130|-||34||. This model is based on the dynamics of a relativistic 
string, representing the colour flux stretched between the initial qq. The string 
produces a linear confinement potential and an area law for matrix elements: 



\M{qq^ h l ---h n )\ 2 oc e 



-bA 



where A is the space-time area swept out (fig. ||). The string breaks up into 
hadrons via qq pair production in its intense colour field. Gluons produced in the 
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12 3 4 

log 10 Mass (GeV/c 2 ) 



Figure 3: Cluster model: mass distribution of qq pairs. 




Figure 4: String model: space-time picture. 
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parton shower give rise to 'kinks' on the string. The model has extra parameters 
for the transverse momentum distribution and heavy particle suppression. It has 
some problems describing baryon production, but less than the cluster model. 



The UCLA model []3~5||-[37| is a variant of the JETSET string model which takes 



the above area law for matrix elements more seriously, using it to determine the 
relative rates of production of different hadron species. This results in heavy 
particle suppression without extra parameters, the mass-squared of a hadron 
being proportional to its space-time area. At present the model still uses extra 
parameters for pt spectra, and again has some problems describing baryon 
production. 



4 Single-particle yields and spectra 

Tables [I] and 2 5 compare predictions of the above models^ with data on Z° decay 
from LEP and SLC Of course, the models have tunable parameters, but the overall 
agreement is encouraging. As stated earlier, the main problems are in the baryon 
sector, especially for HERWIG. 

It is remarkable that most measured yields (except for the 0~ mesons, which have 
special status as Goldstone bosons) lie on the family of curves 



(n) = a{2J + l)e 



-M/T 



where M is the mass and T ~ 100 MeV (fig. |5| f39|). This suggests that mass, rather 
than quantum numbers, is the primary factor in determining production rates. Note 
that, surprisingly, the orbitally-excited J = | baryon A(1520) (not yet included in 



models) is produced almost as much as the unexcited J = | baryon £(1385) 0, [41 



At other energies, model predictions for identified particle yields are in broad 
agreement with e + e~ data (fig. ^| Hj), but statistics are of course poorer. Charged 
particle spectra at low x agree well with the resummed (MLLA) predictions |15, [H], [17 
over a wide energy range, as illustrated in fig. ^ 



In pp — > dijets [fPj the relevant scale is taken to be Q = Mjj sin 9 where Mjj 
is the dijet mass and 9 is the jet cone angle (fig. |8|). Results are then in striking 
agreement with theory and with data from e + e~ annihilation at Q = -^i (fig. [^). 

New SLD data include hadron spectra in light quark (rather than antiquark) 
fragmentation, selected by hemisphere using the SLC beam polarization ]4"5| . One 
sees strong particle/antiparticle differences in the expected directions (fig. [10]), bearing 
in mind the predominance of down-type quarks in Z° decay. 

4 Updated from ref. |§. 

5 Recent ALEPH HERWIG tuning with strangeness suppression 0.8 pSj. 
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Particle 


Multiplicity 


HERWIG 


JETSET 


UCLA 


Expts 






5.9 


7.4 


7.4 




Charged 


20.96(18) 


20.95 


20.95 


20.88 


ADLMO 


7T ± 


17.06(24) 


17.41 


16.95 


17.04 


ADO 


7T° 


9.43(38) 


9.97 


9.59 


9.61 


ADLO 


V 


0.99(4) 


1.02 


1.00 


0.78 


ALO 


,9(770)° 


1.24(10) 


1.18 


1.50 


1.17 


AD 


w(782) 


1.09(9) 


1.17 


1.35 


1.01 


ALO 


r7'(958) 


0.159(26) 


0.097 


0.155 


0.121 


ALO 


fo(980) 


0.155(8) 


0.111 


~0.1 


— 


ADO 


a (980) ± 


0.14(6) 


0.240 


— 


— 


O 


0(1020) 


0.097(7) 


0.104 


0.194 


0.132 


ADO 


f 2 (1270) 


r\ -100/1 A \ 

0.188(14) 


0.186 


~ 0.2 




ADO 


f 2 (1525) 


0.012(6) 


0.021 






D 


K± 


2.26(6) 


2.16 


2.30 


2.24 


ADO 


K° 


2.074(14) 


2.05 


2.07 


2.06 


ADLO 


K*(892) ± 


0.718(44) 


0.670 


1.10 


0.779 


ADO 


K*(892)° 


0.759(32) 


0.676 


1.10 


0.760 


ADO 


K?j(1430) 


0.084(40) 


0.111 






DO 


D± 


0.187(14) 


0.276 


0.174 


0.196 


ADO 


D° 


0.462(26) 


0.506 


0.490 


0.497 


ADO 


D*(2010) ± 


0.181(10) 


0.161 


0.242 


0.227 


ADO 


D± 


0.131(20) 


0.115 


0.129 


0.130 


O 


B* 


0.28(3) 


0.201 


0.260 


0.254 


D 




0.118(24) 


0.013 






D 


J/v 


0.0054(4) 


0.0018 


0.0050 


0.0050 


ADLO 


V>(3685) 


0.0023(5) 


0.0009 


0.0019 


0.0019 


DO 


Xcl 


0.0086(27) 


0.0001 






DL 



Table 1: Meson yields in Z° decay. Experiments: A=Aleph, D=Delphi, L=L3, 
M=Mark II, 0=Opal. Bold: new data this year. Underlined: disagreement with 
data by more than 3a. 
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Particle 


Multiplicity 


HERWIG 


JETSET 


UCLA 


Expts 






5.9 


7.4 


7.4 




P 


1.04(4) 


0.863 


1.19 


1.09 


ADO 


A++ 


0.079(15) 


0.156 


0.189 


0.139 


D 




0.22(6) 


0.156 


0.189 


0.139 





A 


0.399(8) 


0.387 


0.385 


0.382 


ADLO 


A(1520) 


0.0229(25) 








DO 


E± 


0.174(16) 


0.154 


0.140 


0.118 


DO 


S° 


0.074(9) 


0.068 


0.073 


0.074 


ADO 


£*± 


0.0474(44) 


0.111 


0.074 


0.074 


ADO 




0.0265(9) 


0.0493 


0.0271 


0.0220 


ADO 


-(1530)° 


0.0058(10) 


0.0205 


0.0053 


0.0081 


ADO 




0.0012(2) 


0.0056 


0.00072 


0.0011 


ADO 


A c + 


0.078(17) 


0.0123 


0.059 


0.026 


O 



Table 2: Baryon yields in Z° decay. Legend as in table 1. 



LEP1 

Total rates 
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M (GeV/c 2 ) 

Figure 5: Particle yields in Z° decay. 
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Figure 6: Particle yields in e + e annihilation. 
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Figure 7: Low-x fragmentation in e + e annihilation. 



9 



Figure 8: Cone angle in pp — > dijets 




Figure 9: Low-x fragmentation in pp — > dijets. 
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Figure 11: Selection of gluon jets by DELPHI 
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Figure 12: Charged particle spectra in quark and gluon jets. 



5 Quark and gluon jets 



DELPHI [4£| select gluon jets by anti-tagging heavy quark jets in Y and Mercedes 
three-jet events (fig. [□]). As expected, the higher colour charge of the gluon {Ca = 3 
vs. Cp = 4/3) leads to a softer spectrum and higher overall multiplicity (fig. 0). In 
general the relative multiplicities of identified particles are consistent with those of 



all charged, with no clear excess of any species in gluon jets (fig. [13]). In particular 
there is no enhanced 0(1020) or 7] production: 



DELPHI |H: N g (<p)/N q (<p) =0.7 ±0.3 
OPAL [gTj: N g (rj)/N q (r]) = 1.29 ±0.11 

OPAL |48| select gluon jets recoiling against two tagged 6-jets in the same hemi- 
sphere. Monte Carlo studies indicate that such jets should be similar to those emitted 
by a point source of gluon pairs. The qualitative message from the data is again clear 
(fig. [14]): Gluon jets have softer fragmentation than light quark jets, and higher mul- 
tiplicity. The precision of the data is now such that next-to-leading order calculations 
of the relevant coefficient functions, taking into account the experimental selection 
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procedures, are needed to check universality of the extracted gluon fragmentation 
function. 

The ratio of gluon/quark multiplicities at low rapidity (large angle) is close to the 
ratio of colour charges r = Ca/Cf = 2.25, in agreement with local parton-hadron 
duality: 

OPAL: r ch {\y\ < 1) = 1.919 ± 0.047 ± 0.095 

Monte Carlo studies suggest that a better measure of Ca/Cf is obtained by 
selecting low-momentum hadrons with relatively large pt (i.e. low rapidity). This 
gives 

OPAL: r ch (p < 4, 0.8 < p T < 3 GeV) = 2.29 ± 0.09 ± 0.015 

DELPHI |49| have observed scaling violation in quark and gluon jet fragmentation 
separately (fig. |T3|) by studying the dependence on the scale 

k h = E jet sin(0/2) ~ \^fsyi 



where 9 is the angle to the closest jet and y% is the Durham jet resolution |5(J at which 
3 jets are just resolved. This is expected to be the relevant scale when y 3 becomes 
small. One sees clearly that there is more scaling violation in gluon jets (fig. [II]). The 
ratio provides another measure of Ca/Cf- 

DELPHI: r sc . vio i. = 2.23 ± 0.09 ± 0.06 
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gure 14: Momentum fraction and rapidity distributions in quark and gluon jets. 
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Figure 16: Logarithmic gradients of quark and gluon fragmentation. 



A crucial point in the DELPHI analysis is that 3-jet events are not selected using 
a fixed jet resolution y cut , but rather each event is clustered to precisely 3 jets. This 
avoids 'biasing' the gluon jet sample by preventing further jet emission above y cnt . 

The same point is well illustrated in analyses of average multiplicities in 2- and 
3-jet events fTj], [52], |53|]. If N qq -(s) is the 'unbiased' qq multiplicity, then in events 



with precisely 2 jets at resolution y cut there a rapidity plateau of length \n(l/y cut ) 
(see fig. ]17|) and the multiplicity is 

N 2 (s,y cut ) ^ N q q(sy cut ) + \n(l/ y cut )N' q -(sy cut ) 

where N'(s) = sdN/ds. Clustering each event to 3 jets we get this multiplicity with 
2/3 in place of y cut , plus an unbiased gluon jet: 

N 3 (s) ~N 2 (s,y 3 ) + ^N gg (sy 3 ) 

Thus one can extract the unbiased gg multiplicity, plotted in fig. |l8| vs. pj ~ ^/syi 



d4[]. The ratio of gg/qq slopes gives yet another measure of Ca/Cf f5 
r mult = 2.246 ± 0.062(stat.) ± 0.080(sys.) ± 0.095(theo. 



6 Current and target fragmentation in DIS 

HI |55|] and ZEUS have studied the distributions of x p = 2\p\/Q in the current 
and target hemispheres in the Breit frame (fig. [19]). 
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y cut y cut ~ 

Figure 17: Rapidity plateau in two-jet events. 





16 



Tl 
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Figure 20: Contributions to the final state in DIS. 



In the current hemisphere one expects fragmentation of the current jet (C in 
fig. EJ]), similar to half an e + e~ event. In the target hemisphere, the contribution Tl 
is similar to C, T2 gives extra particles with x p < 1, while T3 gives x p ^l, generally 
outside detector acceptance. 

• In the current hemisphere the charged multiplicity is indeed similar to e + e~ 
(fig. p_I] ||56|| ). Differences at low Q 2 are consistent with the expected boson- 
gluon fusion contribution. The distribution of £ = ln(l/x p ) is also similar to 
e + e~, i.e. close to Gaussian with little Bjorken x dependence (fig. p2|) . 

At low Q 2 there is evidence of strong subleading corrections. The distribution 
is skewed towards higher values of £ (smaller x p ), contrary to MLLA predictions 
(fig. E3I). The quantity plotted is 



Skewness = ((£-£) 3 )/((£-0 



2\ 2 



On the other hand, the data lie well below the fixed-order perturbative prediction 
[[57] at low x p and Q 2 (fig. Discrepancies could be due to power-suppressed 
(1/Q 2 ) corrections, of dynamical and/or kinematical origin. The bands in fig. |4] 
correspond to an ad-hoc correction factor 



(0.1 < m cS < 1 GeV). 
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Figure 21: Charged multiplicity in current hemisphere. 



ZEUS 1995 



Q 2 (GeV 2 ) 
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Figure 22: Fragmentation in DIS. Upper data (heavy curve) target region, lower data 
(light curve) current region. 
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ZEUS preliminary 1994-1997 




Figure 23: Skewness in current fragmentation region. 




Figure 24: Scaling violation in DIS fragmentation. 
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Figure 25: Target fragmentation compared with models. 



In the target hemisphere there is also disagreement with MLLA ||56|| , possibly 
due to the T3 contribution "leaking" into the region x p < 1. If anything, 



Monte Carlo models predict too much leakage (fig. ^5]). Little Q 2 dependence 
is evident. 



7 Heavy quark fragmentation 



New data on b — ► B fragmentation from SLD [58], using high-precision vertexing, 
discriminate between parton-shower plus hadronization models (fig. |26|). Note that 
the data have not yet been corrected for detector effects. 

Including more perturbative QCD leads to a reduction in the amount of non- 
perturbative smearing required to fit the data. Non-perturbative effects are conven- 
tionally parametrized by e b in the Peterson function []59|1 

f(z) = - (l - - - t^—) (z = x B /x b ) 

z V z 1 — zj 

Pure Peterson |5_8J: e b = 0.036. 

JETSET (~ LLA QCD) + Peterson @: e 6 = 0.006. 
NLLA QCD + Peterson |6D[: e b = 0.002 (fig. g7|). 
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Figure 26: SLD data on 6 — ► B fragmentation compared with models. 
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Figure 27: SLD []5S| and ALEPH [51] data on b — > B fragmentation, the latter 
compared with NLLA QCD p[. 
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Figure 28: Bose-Einstein correlations with repect to axes along and perpendicular to 
the thrust axis. 



In the universal low-scale a s model, the perturbative prediction is extrapolated 
smoothly to the non-perturbative region, with no Peterson function at all [ pOf| . 



8 Bose-Einstein correlations 



Studies of tt it correlations which distinguish between directions along and perpen- 
dicular to the thrust axis find definite evidence for elongation of the source region 
along that axis (fig. |28| and table |3| |62|, 63], 



This has a good explanation in the 
Lund string model, in terms of the change of the space-time area A in fig. |] when 



identical bosons are interchanged |65 



ALEPH has clear evidence of Fermi-Dirac anticorrelation in AA (S=l). 
A-C in fig. ^9] correspond to different comparison (no-correlation) samples, 
source size appears to decrease with increasing particle mass (table |] |66| 



Plots 
The 



However, some of this effect is kinematic ^9 . 



22 



Expt 


Rl (fm) 


R T (fm) 


Rl/Rt 


DELPHI 
L3 
OPAL 


0.85 ± 0.04 
0.74 ±0.04 
0.935 ± 0.029 


0.53 ± 0.04 
^fi+ - 03 

U - OD -0.06 

0.720 ± 0.045 


1.61 ±0.10 
1.23 ± 0.03±g"f§ 
1.30 ±0.12 ' 



Table 3: Longitudinal and transverse source radii. 
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Particles 


^source (fni) 


7T7T 

KK 
AA 


0.65 ±0.04 ±0.16 
0.48 ± 0.04 ± 0.07 
0.11 ±0.02 ±0.01 



Table 4: Comparison of source radii. 




Figure 30: Hadronizations volumes in WW decay. 



9 WW fragmentation 

In e + e~ — ► WW, we would expect correlations between W hadronic decays due to 
overlap of hadronization volumes. This occurs mainly in the central region, and is 
orient at ion- dependent (fig. ^). These reconnection effects have been searched for in 
single-particle distributions. One would expect discrepancies between the distribu- 
tions in semi-leptonic and fully hadronic decays, especially at low momenta. There 



is no firm evidence yet for such effects in the x p distribution (fig. |3_T] [70]). However, 



DELPHI [Q report a possible small (~ 2a?) effect in the distribution of px relative 



to the thrust axis (fig. [32]) 

Bose-Einstein correlations between hadrons from different Ws are also being 
looked for. They would lead to an increase in the correlation function for WW 
relative to that for a single W. There is no sign of any increase at present (fig. ^ 

OH- 



IO Summary 

• Detailed fragmentation studies need more theoretical input in the form of coef- 
ficient functions that take account of selection procedures, especially for gluon 
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ALEPH PRELIMINARY 189GeV 




Figure 31: Distribution of momentum fraction in WW decay. 
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;ure 32: Distribution of transverse momentum in WW decay. 
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A = 0.23 ±0.03 

a = 4.26 ± 0.43 GeV"' 
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Figure 33: Bose-Einstein correlations in WW decay. 
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jets in e + e final states. 

• Hadronization studies suggest that particle masses, rather than quantum num- 
bers, are the dominant factor in suppressing heavy particle production. Baryon 
production is not yet well described by any model. 

• Quark and gluon jets have the expected differences and these can be used to 
measure the ratio of colour factors Ca/Cf- There is no strong evidence yet for 
different particle content in gluon jets. 

• Fragmentation in DIS shows disagreements with perturbative predictions. It is 
not yet clear whether these are due to higher-order or non-perturbative effects. 

• New precise b quark fragmentation data test models and suggest that pertur- 
bative effects dominate. 

• Bose-Einstein (Fermi-Dirac) correlations show elongation of the source along 
the jet axis and source shrinkage with increasing mass. 

• WW fragmentation still shows no firm evidence for correlation between the 
decay products of the two Ws. 
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Discussion 



Charles Buchanan (UCLA): To elaborate on the UCLA approach to hadroniza- 
tion: We find that the idea of a spacetime area law (as suggested by strong QCD) 
works very well as an organizing principle in the stage of the soft strong-coupled 
hadronization part of the process — that is, it gracefully predicts "easy" data such as 
light-quark meson production rates and distributions in e + e~ with few parameters 
and forms an attractive basis for studying the relation with the perturbative stage 
and more complicated phenomena such as baryon formation and pr effects. To study 
these latter in detail, we have joined BaBar where we will use the 10 8 high quality 
events to be collected in the next 2-3 years to study for baryon-meson-antibaryon 
3-body corrections, pr correlations, etc. We invite physicists interested in the area to 
contact us (at buchanan@physics.ucla.edu). 

George Hou (National Taiwan University): Regarding particle content of gluon 
jets, is there any result on the r( content? This particle is more naturally associated 
with gluons than rj or 0. 

Webber: As far as I know, there are no results available yet on the rj' content of 
gluon jets. 

Michael Peskin (SLAC): You have shown that, when quark and gluon jets are 
carefully selected, their average properties are clearly distinguished. Of course, what 
one really wants is a variable which allows one to separate quark and gluon jets (if 
only statistically) in the Tevatron or LHC environment. What is the best choice for 
this purpose? 

Webber: This is difficult to do because, although the average properties are differ- 
ent, the fluctuations are large. 

Tom Ferbel (University of Rochester): Actually, DO has used the differences 
between quark and gluon jets very effectively, in a statistical manner (with neural 
networks) to improve the signal to background in the analysis of it production in the 
all-jets channel. 



31 



